Sketching Non-linear Systems 



In session on Phase Portraits, we described how to sketch the trajecto- 
ries of a linear system 

x' = ax + by , , ^ ^ 

, , / a,b,c,a constants. 

y = cx + ay 

We now return to the general (i.e., non-linear) 2x2 autonomous system 
discussed at the beginning of this chapter, in sections 1 and 2: 

x' = f{x,y) _ 
y' = g{x,y) ' 

it is represented geometrically as a vector field, and its trajectories — the 
solution curves — are the curves which at each point have the direction 
prescribed by the vector field. Our goal is to see how one can get informa- 
tion about the trajectories of (1), without determining them analytically or 
using a computer to plot them numerically. 

Linearizing at the origin. To illustrate the general idea, let's suppose that 
(0, 0) is a critical point of the system (1), i.e., 

/(0,0)=0, ^(0,0) =0, (2) 

Then if / and g are sufficiently differentiable, we can approximate them 
near (0, 0) (the approximation will have no constant term by (2)): 

f{x, y) = uix + hiy + higher order terms in x and y 
g{x, y) = a2X + bzy + higher order terms in x and y. 

If (x, y) is close to (0, 0), then x and y will be small and we can neglect the 
higher order terms. Then the non-linear system (2) is approximated near 
(0, 0) by a linear system, the linearization of (2) at (0,0): 

x' = aix + biy 
y' = + b2y ' 

and near (0,0), the solutions of (1) — about which we know nothing — will 
be like the solutions to (4), about which we know a great deal from our 
work in the previous sessions. 

Example 1. Linearize the system < ^ , ^^^^^ at the critical 

^ \y' =x{l+yY 

point (0,0). 
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Solution. We have < ^ 1''^ ^ so the linearization is 

x' = y 
y' = X ' 

Linearising at a general point More generally, suppose now the critical 
point of (1) is (xo,yo)/ so that 

/(^o,yo)=0, ^(xo,yo)=0. 

One way this can be handled is to make the change of variable 

x-y = x-xq, yi = y-yo; (4) 

in the Xii/i-coordinate system, the critical point is (0, 0), and we can proceed 
as before. 

f ^ X — x^ — 2x1/ 
Example 2. Linearize < , 9 q at its critical points on 

{ y =y-y -jxy 

the X-axis. 

Solution. When y = 0, the functions on the right are zero when x = and 
X = 1, so the critical points on the x-axis are (0,0) and (1,0). 

The linearization at (0,0) is x' — x, y' = y. 

To find the linearization at (1,0) we change of variable as in (4): Xi = 
X — 1, yi = y ; substituting for x and y in the system and keeping just 
the linear terms on the right gives us as the linearization: 

x[ = (xi + l)-(xi + l)2-2(xi + l)yi -xi-2yi 

y'l = yi - y? - 1(^1 + i)yi ~ -|yi • 

Linearization using the Jacobian matrix 

Though the above techniques are usable if the right sides are very sim- 
ple, it is generally faster to find the linearization by using the Jacobian ma- 
trix, especially if there are several critical points, or the functions on the 
right are not simple polynomials. We derive the procedure. 

We need to approximate / and g near (xq, yo)- While this can sometimes 
be done by changing variable, a more basic method is to use the main ap- 
proximation theorem of multivariable calculus. For this we use the notation 

Ax = x-xo, Ay = y-yo. A/ = /(x,y) - /(xo,yo) (5) 
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and we have then the basic approximation formula 

since by hypothesis /(xq, yo) = 0. We now make the change of variables (4) 

Xi = X-Xq = ^X, y-^ = y-yo = Ay, 

and use (6) to approximate / and g by their linearizations at (xo,i/o)- The 
result is that in the neighborhood of the critical point (xo,yo), the lineariza- 
tion of the system (1) is 



(7) 



= Axi, where xi = ( ^ and A = f 



In matrix notation, the linearization is therefore 

^0 and 

(8) 

the matrix A is the Jacobian matrix, evaluated at the critical point (xo,yo)- 

General procedure for sketching the trajectories of non-linear systems. 

We can now outline how to sketch in a qualitative way the solution 
curves of a 2 x 2 non-linear autonomous system, 

y' = 8{x>y)- 

1. Find all the critical points (i.e., the constant solutions), by solving the 
system of simultaneous equations 

/(x,y)=0 
g{x,y) = 0. 

2. For each critical point (xq, yo), find the matrix A of the linearized system 
at that point, by evaluating the Jacobian matrix at (xo,yo): 



Sketching Non-linear Systems 



OCW 18.03SC 



(Alternatively, make the change of variables Xi = x — 
^0/ yi = y ~ yO' ^rid drop all terms having order higher 
than one; then A is the matrix of coefficients for the linear 
terms.) 

3. Find the geometric type and stability of the linearized system at the 
critical point point {xo,yo), by carrying out the analysis in sections 4 and 5. 

si The subsequent steps require that the eigenvalues be non- 
zero, real, and distinct, or complex, with a non-zero real 
part. The remaining cases: eigenvalues which are zero, re- 
peated, or pure imaginary are classified as borderline, and 
the subsequent steps don't apply, or have limited applica- 
tion. See the next section. 

4. According to the above, the acceptable geometric types are a saddle, 
node (not a star or a defective node, however), and a spiral. Assuming 
that this is what you have, for each critical point determine enough addi- 
tional information (eigenvectors, direction of motion) to allow a sketch of 
the trajectories near the critical point. 

5. In the xy-plane, mark the critical points. Around each, sketch the trajec- 
tories in its immediate neighborhood, as determined in the previous step, 
including the direction of motion. 

6. Finally, sketch in some other trajectories to fill out the picture, mak- 
ing them compatible with the behavior of the trajectories you have already 
sketched near the critical points. Mark with an arrowhead the direction of 
motion on each trajectory. 

If you have made a mistake in analyzing any of the critical 
points, it will often show up here — it will turn out to be im- 
possible to draw in any plausible trajectories that complete 
the picture. 

Remarks about the steps. 

1. In the homework problems, the simultaneous equations whose solutions 
are the critical points will be reasonably easy to solve. In the real world, 
they may not be; a simultaneous-equation solver will have to be used (the 
standard programs — MatLab, Maple, Mathematica, Macsyma — all have 
them, but they are not always effective.) 

2. If there are several critical points, one almost always uses the Jacobian 
matrix; if there is only one, use your judgment. 

3. This method of analyzing non-linear systems rests on the assumption 
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that in the neighborhood of a critical point, the non-linear system will look 
like its linearization at that point. For the borderline cases this may not be 
so — that is why they are rejected. The next two notes explain this more 
fully 

If one or more of the critical points turn out to be borderline cases, one 
usually resorts to numerical computation on the non-linear system. Occa- 
sionally one can use the reduction to a first order equation: 

dx fix,y) 
to get information about the system. 
Example 3. Sketch some trajectories of the system 

x' = —x + xy 
y' = -2y + xy ' 

Solution. We first find the critical points, by solving 

—x + xy =x{—l+y)=0 
-2y + xy =y(-2 + x) =0 ' 

From the first equation, either a: = or y = 1. From the second equation, 
x = Q ^y = 0; y = l =^x = l; critical points : (0,0), (2,1). 
To linearize at the critical points, we compute the Jacobian matrices 

i-y y -2 + xJ' ho,o)-[ -2 J h2,i)-[i 

Next we analyze the geometric type and stability of each critical point: 

(0,0): 

eigenvalues: ti = — 1, = —2 sink node 

eigenvectors: ^i~(^o)' 

By the node-sketching principle, trajectories follow eii near the origin, 
are parallel to a.2 away from the origin. 

(2,1): 

eigenvalues: ti = \/2, ^2 = — unstable saddle 

v/2 \ . ^ / -v/2 



eigenvectors: oi\= [ j ) ' '"''^ \ \ 
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Draw in these eigenvectors at the respective points (0, 0) and (2, 1 ), with 
arrowhead indicating direction of motion (into the critical point if 1 < 0, 
away from critical point if 1 > 0.) Draw in some nearby trajectories. 

Then guess at some other trajectories compatible with these. See the 
figure for one attempt at this. Further information could be gotten by con- 
sidering the associated first-order ODE in x and y. 




Example. Sketch the phase portrait of the following system. 



14x 
161/ 



xy 
xy 



Critical points: 

x(u-'^x-y\ = {) 

y ^y-^j =0 



or 14 X - 1/ = 

2 



y = or 16 



x = Q =^ y = or y = 32. 

y = =^ X = or X = 28. 

x^O,y ^0 => X = 12, y = 8. 

=> all critical points: (0,0), (0,32), (28,0), (12,8). 



0. 
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14 — X — y —X 

V -y 16-1/-X ^ 

Looking at each of the critical points in turn: 

7(0, 0) = ^ 16 ) ■ ^^S^^^^^^^^ 14, 16; eigenvectors ^ g ) ' 
=^ source node (see 'Source node' picture below). 

32 16 ) " ^'^S^^^^'^"^^^ "1^' '1^' eigenvectors ( 
^ sink node (see Sink node 1' picture below). 

7(28,0) = 

^ sink node (see Sink node 2' picture below). 

7(12,8) = (:| 

eigenvalues — 5 ± \/97 ^ —15,5; 
eigenvectors I ^ 1,1 

^ saddle (see 'Saddle' picture below). 

Rough sketch of system: 
First we sketch each of the critical points. 

V 



-12 
-4 



11 
8 






Tr / 







Source node 

V 



Sink node 1 
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Sink node 1' picture below). 

-14 -28 \ . , ..... f 
Q 12 j ' ^^S^^'^^^^^^ "14, -12; eigenvectors I 

CI 1 3 „ • _j 1 1 \ 



1 
16 

1 ^ 

y 




Sink node 2 



Saddle 
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Computer plot - phase plane portrait. 
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